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Liouville ergodicity of linear multi-particle hamiltonian 
system with one marked particle velocity flips 

Lykov A. A., Malyshev V. A.* 


Abstract 

We consider multi-particle systems with linear deterministic hamiltonian dynamics. 
Besides Liouville measure it has continuum of invariant tori and thus continuum of invari¬ 
ant measures. But if one specified particle is subjected to a simple linear deterministic 
transformation (velocity flip) in random time moments, we prove convergence to Liouville 
measure for any initial state. For the proof it appeared necessary to study non-linear 
transformations on the energy surface. 
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1 Introduction 


Ergodicity problem for hamiltonian multi-particle systems produced many deep results. First 
of all, many examples of non-ergodic systems appeared - linear, non-linear with additional 
integrals and close to them (KAM theory). One could expect then that for generic hamiltonians 
one also has non-ergodic behaviour. However, as far as we know, this is still an open difficult 
problem with many partial results, see [201119], and, after the century history of the ergodicity 
hypothesis it is reasonable to look for simpler alternative approaches to it. 

Namely, one could assume that any physical system has always some contact with external 
world. Such contact can be of quite various extent: 1) all particles can have contact with 
external world and stochastic behaviour (for example, with dynamics of Glauber type), 2) only 
particles on the boundary, etc. But then it is quite natural to ask - what is the minimal 
contact which dehnitely provides ergodic behavior. Possible reformulation of the ergodicity 
hypothesis could be the following; for generic system even the minimalistic contact produces 
ergodic behaviour. 

Here we consider an example of such minimalistic contact which consists, hrst of all, in 
that we allow some contact with external world for only one (marked) particle. In our earlier 
papers mu this particle was subjected to some random force, that garantied convergence 
to Gibbs equilibrium. Here we assume even less randomness. Namely, the marked particle is 
subjected to a simple deterministic transformation (velocity flip, very popular in other problems 
[a m El El El m) but in discrete random time moments 

0 < ti < ... < tjn < ■ ■ ■ ( 1 ) 

Liouville measure is evidently invariant w.r.t. such dynamics on the energy surface. We prove 
then that ergodicity in the stronger form holds - for any initial state we have convergence to 
Liouville measure on the energy surface. 

It is very interesting that it works even for the systems having the worst possible non- 
ergodic behaviour - linear systems. The only price we pay is that ergodicity holds not for any 
linear system but for almost any - this is purely algebraic phenomenon which was discussed 
in our earlier papers and cannot be avoided. It seems reasonable that the same result holds 
for non-linear systems as well (following the common belief that the latter have better mixing 
properties than linear systems). 

The paper is naturally subdivided in two parts. The hrst part uses no probability at all but 
only elaborate non-linear analysis to prove that the trajectory visits all invariant tori and even 
any point - we had to use coordinates on the energy surface where velocity hips are strongly 
non-linear. Second part, on the contrary, essentially uses non-trivial parts of Markov processes 
theory on continuous state space. 

Now we come to rigorous dehnitions. 


2 Model and Main Results 

Hamiltonian dynamics We consider the linear space 

L = = {^l; = : q= (gi,...,gjv)^, P = (pi,..., pjv)^, qi,Pi e M}, 


where T denotes transposition (thus -0 is a column-vector). It can be presented as the direct 
sum L = 0 of two orthogonal (coordinate and momenta) spaces of dimension N with 


the standard scalar product in M' 


2N 


N 

= (g,g')2 + {p,p')2 = ^{qiq'i + PiP'i) 

i=l 

We consider quadratic haniiltonian 


^^2 1 

k=l 

where the matrix 1/ > 0 acting in is assumed to be real and positive dehnite (thus the 
particles cannot escape to inhnity). This dehnes hamiltonian system of linear ODE with k = 

N 

qk = Pk,Pk = - Vkiqi, (3) 

i=i 

For any h > 0 dehne the constant energy surface 

Mh = {^l’eL: = h} 

Then AAh is a smooth manifold (ellipsoid) in L of codimension 1. 

Allowed hamiltonians Dehne the mixing subspace 

L_=L_(r) = {f O eL:q,pelv} 

where ly = lv,i is the subspace of generated by the vectors V^ei, k = 0,1, 2..where 
Cl,C at is the standard basis in R^. 

Let V be the set of all positive-dehnite {N x A)-matrices, and let V'*' C V be the subset 
of matrices for which 


L-{V) = L (4) 

The set of D G V such that their eigenvalues, denoted by ojf,... , 00 %, are independent over the 
held of rational numbers, is denoted by Vind- 

Lemma 1 The set is open and everywhere dense (assuming topology of R 2 ) in V, 

and the set V+ 0 ^ind is dense both in V+and in V. 

See more in section 5.3. 


Piecewise deterministic process Assume that at time moments ([I]) the following deter¬ 
ministic transformation / ; L —)■ L occurs; all qkiPk are left unchanged, except for pi, the sign 
of which becomes inverted 


Plitm - 0) -^Pl(tm) = -Plitm - 0),m > 1 

For example, one can consider L as the phase space for N identical point particles in R, with 
mass m = 1, and real numbers are their coordinates and velocities (momenta). Then 



this transformation can be interpreted as the elastic collision of the particle 1 with a wall. 
Alternatively, taking dN instead of N, one can imagine N particles in R‘^ where only one 
velocity component of particle 1 is flipped. Reflections w.r.t. any hyperplane in R‘^ could be 
considered quite similarly. 

In-between these moments the system evolves via hamiltonian dynamics Q. 

With {2N X 2iV)-matrix 


A 


0 E 
-V 0 


the system ([2]) can be rewritten as 


ij) = All;. 


( 5 ) 


and the solution '0(t) of ([2]) with initial vector '^(0) will be 


= e^^7p{0). 


For given sequence ([I]) the dynamics of our process is defined for tm < t^+i as 


where ti = ti,T2 = t2 — ti, ..., Tm = tm — tm-ii •••• For any f ^ 0 define linear maps L ^ L 

^ L 

It is clear that AA.h is invariant w.r.t. J{t) for any h > 0 and f > 0. For any 'ip & L and any 
integer m ^ 1 define the set of states 


Jmi.'pj) = ■ ■ ■ J{.n)'lp : 0 < Ti, C Mh 

which the system can visit at the m-th flip. 

Theorem 1 (covering theorem) Assume that V G V”*" fl then there exists m ^ 1 such 
that for any xp & L we have 

PTmidP^ ~ -Al/i 


We introduce randomness by the following assumption 
AO); positive random variables 


Ti fi, T2 ^2 ^1; ■■•5 T~n i'n 15 ••• 

are assumed to be independent, identically distributed with measure Pr = p{s)ds where the 
density p (w.r.t. Lebesgue measure ds) is positive everywhere on i?+, and moreover the first 
moment Eti < oo. 

If, for example, r* have exponential distribution with the density Aexp(—Ar), A > 0, then it 
defines Markov process -^(f) with right continuous deterministic trajectories and random jumps. 
Such processes are often called piecewise deterministic Markov processes, see for example ig. 
At the same time, this can be considered as an example from random perturbation theory, see 
ra where the problem of invariant measures is studied. 

Let TT be Liouville measure on the energy surface, defined by the surface form da divided by 
|ViL|. It is well known that tt is invariant w.r.t. hamiltonian dynamics and also w.r.t. velocity 
flips. 

Theorem 2 (convergence theorem) Assume that V G H Then under assumption 

AO) for any initial '0(0) and any hounded measurable real function f on Adh we have a.s. 

MjiT) P I nm)dt ^ ir(/) [ fdi, 

Jo JMh 

We call this property Liouville ergodicity. 


3 Deterministic part - proof of covering theorem 

Plan of the proofs In all assertions below we always assume that V G H Theorem 
[2] will follow from theorem [H Theorem [T] will follow from the following weaker results. 

Theorem 3 (closure theorem) There exists m ^ 1 such that for all & L we have 

Jmifj) = Mh 

The covering and closure theorems have simpler local analogue. And moreover it shows 
that the dimension of ffmi'ip) grows as m for m = 1, 2,..., 2N — 1. For exact formulation we 
need some dehnitions. For any ri,..., Tm > 0 denote J(ri,..., r^) = Jijm) ■ ■ ■ Jiji)- It can be 
considered as the mapping from L to M.h for fixed (ri,..., Tm), but also J(ri,..., Tm)'ip can be 
considered as the map 

Mh, f = (n, 

from the m-dimensional orthant Qm = {(r'l,..., r^) : Xj > 0, i = 1,..., m} = to Mh, for 
fixed Ip e Mh- 

Theorem 4 (local covering theorem) For any point ip G Mh and any k = 1,, 2N — 1 

the dimension of J^flk is equal to k. Moreover, there exists open subset U <zVtk such that the 
mapping : U ^ Mh is a smooth embedding 

3.1 Key definitions, notation and some intuition 

The sequence ([T]) completely defines the trajectory, or the path. Proof of convergence and 
covering theorems could be obtained, in some sense, by “summation” over all possible paths. 
For this we have to use various coordinate systems on the energy surface. To get some intuition 
it is useful to see how it works for easier cases N = 1,2. 

Action-angle coordinates Let ui,..., uat be the orthonormal eigenvectors of V in and 
let ujf ,..., Lv^ be the corresponding eigenvalues of V. For any k = 1,..., N define 2A-vectors 
in L: 

Qk = (vk,0)'^, Pk = (0,Vk)'^ (6) 

The coordinates of ■0 G L in the basis Qi, Pi,..., Qjv, Pn we denote by {qk,Pk)'^, that is 

N N 

V' = X] QkQk + ^ PkPk (7) 

k=l k=l 

They correspond to coordinates and momenta of the independent «quasiparticles», that is one- 

r'i 

dimensional oscillators, having energies where 

rlW =pI+ ojIqI = (V’, Pk)l + ujli'ip, Qk)l = {p, Pk)l + Qk)l, k = l,...,N 

are the action coordinates of the point tp = {q,pY. We agree that Vkigp) = ■\/r^('^) ^ 0. It is 
easy to see that are integrals of the hamiltonian dynamics, that is for any f ^ 0 

( 8 ) 

The angle variables then are the angles for these oscillators. 

The following assertions easily follow from the known facts, see for example [T], pp. 103, 
272. But the proof is very elementary and we give it for the reader’s convenince in Appendix. 




Lemma 2 Consider hamiltonian dynamics then 


1. For any ri ^ 0,..., rjv ^ 0 the set 

T(ri, ...,rN) = {^jJEL: rk{ip) = rk,k = 1,A^} 

is invariant and diffeomorphic to torus of dimension N — n, where n equals the number 
of zeros among ri,..., Ttv- 

2. For any point E L the closure of its orbit coincides with the torus defined by it, that is 

t ^ 0} = T{ri{'ip ),..., VNi'ip)), 

3. Thus the torus defines the vector f = ri, r 2 , • • •, tat. Vice-versa, any such vector with 
non-negative coordinates uniquely defines the torus. This torus lies on the energy surface 
M, zjf 

N 

'^rl = 2h 

k=l 


For convenience we put h = };, and denote the set of all invariant tori on Ad = Adi by 

^ 2 

N 

T = {f = (ri,..., r^v) : T-i ^ 0, ^ = 1} 

k=l 


Case N = 1 Here uf can be arbitrary, and we take h = = 1, let S' be the corresponding 

circle. Then the particle moves along S with constant angle velocity in the clock-wise direction. 
Covering theorem is evident in this case (with m = 1). However, for = 1 we have much 
stronger statement. 

For any denote T{'ijj,ijj') the minimal t > 0 such that ip' = Then Tq = T{ip,ip) 

is the first return time, or the time of complete rotation around the circle, and for any ip 

T{iP,-iP) = ]^T, 

Proposition 1 For any ip, ip' and any t > Tq there exists 0 < = ti{ip,ip',t) < Tq such that 

= iP' 


For the proofs see Appendix. 

(r, p)-coordinates and (r, p)-coordinates Note that I can be written as 

I = E-2Pi, 

where E is the identity matrix and Pi is the orthogonal projector on the vector gi = (0, ei)^ E L. 
The expansion 

N 

ei = ^ /3kVk 

k=l 



of the vector ei G defines the numbers fik- The pk and pk coordinates of the vector p in 
are related by the formulas 

N N 

p== y^jkVk 

i=l k=l 

N N 

Pk = '^Pi{ei,Vk)2, Pk = '^Pi{vi,ek)2 (9) 

i=l i=l 

Note also that I acts only on p^-coordinates and p^-coordinates, but not on and qk- Thus 
the following notation is justihed 


N 


Ip^p- 2(p, eijjci = p't = Pk- 2piA 


( 10 ) 


k=l 


where 


N 


Pi = 


(P,ei)2 = ^Pk/3k- 


k=l 


By ([9]) the velocity flip changes the action variables as follows 

= (M Vk)l + Vk)2 = ip'kf + ^k)l =Pl + ^plf^l - 4pfcPi/?fc + Vk)l = 

= rlpp) + 4p2/?2 - ApkPil3k = rl{il)) + Apll3l - Api/3k{p, Vk)2 
The following formula dehnes the time evolution of the tori, in terms of initial action variable 
and momenta at time t — 0. Namely, for any k 

= rl{Rli{t - 0)) = rliij) + 4pl{t - 0)/?| - 4pi(t - 0)/3k{p{t - 0),Vk)2, (H) 

where 'ip{t — 0) = 6*^-0 = {q{t — 0),p(f — 0))"^ G L, p{t — 0) = (pi(f — 0),... ,PAr(t — 0))"^. We 
can rewrite OB as 

= T(f(?/>),p(f-0)), (12) 

where 

T = (Ti,...,Tjv),T^(f,p) = rl + Apll3l - Api(3k{p,Vk)2 

has the domain of dehnition 

I1(T) = {(r,p) : f G T,p G Tp{f)} C T x 

where Tp{f) is the projection of the torus T(f) onto the space of momenta, that is the set of 
all p G such that for some q G we have {q,p)'^ G T(f). For any ip = V’(O) = {q,p)'^, q = 
g(0),p = p(0) consider the following set of tori 


(13) 


Tiipj) = {r{J{t)'ip) : f ^ 0} = {A/{r{'ijj),p{t - 0)) : t > 0,p(0) = p} C T. 

From continuity of T it follows that Tj ('?/’) depends only on the invariant torus containing 
Ip, that is if fpip) = fpip') for two points ip^ip' G L, then 

tm = r,{ip') 

Then, the following notation is correct 

T,{iP) = T,{fW) 

It is useful to note that by continuity of T for any f G T we have 

7j(f) = {T(f,p) : p G Tp{f)} 


(14) 




Case N = 2 Here we will only give intuitive arguments why for any initial state 'i/j, in finite 
number of jumps, we can enter any torus. From flTTl) one can get the following formula for the 
evolution of the action variables 

rKJlt)^}) - ri(V’) = 4/3 i/?2(Pi(^ - 0)Pi+p2{t - 0)/32){p2{t - 0)/?i - pi{t - Oj/^a), 
and also 

Denote the right-hand side of the first formula by D{t). It appears that for some t' the set 
{D{t) ■. 0 < t < t'} contains an open interval, and moreover its length has lower bound <5 > 0 
uniform in ip. Thus the measure of the set of visited tori, after each application of I enlarges by 
additive constant. It follows that all tori may be visited for finite number of I transformations. 
To see this for general N is more difficult. 


3.2 Proof of closure theorem 


3.2.1 Contraction property 


Define the function p in T 


N 

p(f,f') = ^\rl 

k=l 



where f = (ri,... ,rAr), f' = (rj,... Denote f* = (|/5i|,..., I/SatI) G T, thus it is the 

invariant torus, containing the point gi = \/2(0,ei)^ G L. Note that all j3k are nonzero, 
that follows from the assumption V G V"*", see section 7.3. The point gi corresponds to the 
configuration of particles where all particle have zero velocity and zero coordinates, except for 
the particle 1. 


Theorem 5 (contraction theorem) For any r = (ri,...,r7v) G T we have the following 
contraction bound 

p(ri(f),f*) ^ (1 - c(f))p(f,f*), 
where the constant c(f) is given by 


1 r-k . Vk 

= - ri ^ = niax ■— - mm — 

max{l,D^(r)} k=i,...,N\f3k\ k=i,...,N \f3k\ 

Proof. We will find point p G Tp{f) such that 

p(T(f,p),f*) = (1 — c(f))p(f, f*). 


(15) 


Note first that in coordinates p the set Tp{f) is the iV-dimensional cube with sides (2ri,..., 2rAr), 
this follows from the oscillator representation. In other words, the point p = (pi,... ,Pn)'^ G 
TJf) iff for any k = 1,..., N 

\Pk\ < Tk, (16) 


For k = 1,..., N denote 


Ik = 


Tk 


and denote the minimal of them by and maximal by y^r. We will need the following functions 







of a; G M, where the constant c is dehned in the formula ffTSj) . As c ^ 1, then for all x 


+ c(l — > 0 

Dehne the point p G with coordinates , pj^) in the basis ui,..., uat by the formula 

« Pk-^l ,+ t ^ 

Pk = y/3k - c , y = ri'yn)- 

and hnd the value of '^{f,p). As 

N N N 

Pi = {p, 61)2 = Y 1 = yJ 2 f^k-^ Y^i^k - rl) = y- 

k=l k=l y k=l 


then for any fc = 1,..., A^ we have 

^1=4 + 4pp| - 4pi&pl = rj + iy^Pl - iyPt (yh - 
Then we have the distance 


rl + c{(3 


2 

k 



N 

k=l 


Lemma 3 p' G Tp{r) that is the inequalities [Wi) hold. 

Proof consists of simple calculations and is given in Appendix. 

From the contraction theorem only, one cannot prove the closure theorem because c = c(f) 
depends on f. That is why we need one more assertion. For any f = (ri,..., tat) G T put 

A(f) = min ^, B{f)= max ^, A{f) = (f) - (f). 

The function A(f) defined on T also looks like a distance to the point f*. Let ns explain this 
in more detail. In particular, A(f) = 0 iff = \f3k\ for all fc = 1,..., A. Moreover, as 

N 

k=l 

then, if A(f) is small then the numbers are close to |/3fc|. Even more, it gives an upper bound 
for the distance, more exactly for any f G T 

p(f,A) = ^ \rl-l3l\ = - 1| ^ l^k{B\f) - A\f)) = '^l3lA{f) < A(f) (18) 

k k ^ k k 

The following result shows that each velocity flip at appropriate moments makes point pi closer 
by 1. 

Corollary 1 For any point f G T there exists p G Tp{f) such that 


A(T(r,p)) ^ max{A(r) — 1,0} 






Proof. From the proof of contraction theorem it follows that there exists p G Tp{r) such 


that 


'^lir.p) = rl + c{r){l3l-rl). 

Denote f' = (r'^,... ,r^) = Without loss of generality one can choose the indices so 


that 


Then 




1/9: 


1/9. 


N\ 




W Pi 


r^N n 


from where we get 
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A(f') ^ (1 — c(f))A(f). 


(19) 


c r = 


max{l, {B{r) — A(r))2} 
Note that the following inequality holds 

{B - Af ^B^-A^ = A. 

and then 


c r = 




max{l, {B{r) — A(r)y} max{l, B‘^{r) — A^{r)} max{l, A(r} 
From this inequality and the bound (IT^ finally we have 

1 


A(r') ^ 1 - 


max{l, A(r)} 


A(r) = max{A(r) — 1, 0}. 


3.2.2 Closure theorem 

Before proving it we have to prove even weaker assertion. 

Lemma 4 There exists integer m ^ 1 such that for any ijj ^ Ai 

gi £ 

In other words there exists m such that for any ip E M. and any e > 0 one can find moments 
0 < fi < ... < so that 

||J(r^)... J(ri)V’ - gi \\2 < e, 

By continuity of maps T and A, and by Corollary [T] it follows that for any point if E L and 
any e > 0 one can find time moment f ^ 0 such that 

A{f{J{t)ijj)) = A(T(f('0),p(t — 0))) < max{A(f('0)) — 1,0} + e. 

It follows that there exists m ^ 1 such that for any e > 0 there exist time moments ti,... ,tm ^ ^ 
such that 


A(r(V’')) ^ e, fj' = J{Tm)...JiTi)7p. 












and moreover one can take m = [A(r('^))] + 1. But as for any r G T the following inequality 
holds 


A(r) ^ max 

k=l,...,N (3l 

m can be chosen uniformly in 'ijj. By formula (IT^ we get 

p(f('0'),f*) ^ e. 


Alternatively it is evident that there exists constant c > 0 such that for any e' > 0 there is 
f ^ 0 such that 

- 9i\\l ^ cp{r{ilj'),f*) + e'. 

As e and e' are arbitrary and by = J(0)J(f) we get the proof. T 
Now we prove the closure theorem. Define another norm on L 

WiPWh = 

Let us fix two points & M. and show that tfj' G for some m ^ 1 not depending on 

By lemma m there exist > 0 and r(, ...,r^ > 0 such that 

- J{t[, . . . < 6 , 

It is clear that the transform J{t[, ... ,r^) is invertible and conserves the norm || \\h, that is 
why 

||J"^(r{,...,r^)J(ri,...,r^)V’-V’llrr < e- 

One can hnd m' ^ 2m and r",..., r^, > 0 such that 

, ■ ■ ■ < 2e, (20) 


In fact 

( 21 ) 

As conserves the norm 11.||ir and / = J(0), h is sufficient to show that for any point %Ij E M., 
any e > 0 and f ^ 0 there exists s = s('^, f, e) such that 


But as the closure of the orbit of any point is the torus, it is Poincare recurrence theorem. 
From (|20|) and equivalence of the norms || ||2 and || ||h the closure theorem follows. 


3.3 Proof of local covering theorem 

Fix xjj. It is sufficient to show that for some point r G Dfc the rank of dJ^{f) equals k. The 
columns of the Jacobian matrix are the vectors 

Okir) = ..., Tk)'ip = /e"''=^AJ(ri,..., Tk-i)'ip 

CLTk 




Denote dim(wi, W 2 ,...) the dimension of the linear span of vectors wi,W 2 ,---- As the map¬ 
ping is non-degenerate, then for any k = 1,2,... one should derive that there is f = 

(ti, ..., Tfc) G Dfc such that 

dim{6'i(r),... ,6'fc(r)) = dim(AJ(fc)'i/', ■ ■ ■, -T^J{rc)'4>) = k, 

CLTi UTfe-l 

where Tc = (ri,..., Tk-i) G Dfc-i- We shall do it by induction in k. Remind that for any 'll) E L 

Ilf) ='ll)- 2 (' 0 , 5 ( 1 ) 25-1. 

In fact, for k = 1 this is trivial. 

Induction step. Assuming that we have proved the assertion for k < 2N — 2, we shall prove 
it for k + 1. 

For f = (ri,..., Tk) G Vtk denote fc = (ri,..., Tk-i) G Dfc_i. Then 

AJ{f)'ip = J(ri,..., Tk-i)'il) = J(ri,..., Tk-i)'ip + ao{f)Agi = 

= (AJ(ri,..., Tk-i)'il) + ao(r)Ae“^'=^ 5 (i) = {AJ{fc)'ilJ + ao{f)Ae~'">^^gi) , 

where ao(r) = —2(e'^'="^J(fc)'0,fi'i) 2 - Similarly we get equalities for the derivatives for i = 

-^J{f)'il) = .. ..Tk-ilip = ... ,rfc-i)V' + ai{r)gi = 

CLTi CLTi CLTi 

= , ai{f) = -2{e^*‘^-^J{fc)'ip,gi)2- 

Derivative in Tk has the following expansion 

= /Ae"'''^J(ri,... = Ae"’^^J{Ti,...,Tk-i)'il) + ak{f)gi = 

CLTk 

= 6^'="^ {AJ{fc)i> + ak{f)e~'^^^gi) , afc(f) = -2(Ae"'''^J(fc)V’, 91 ) 2 - 
Then by nondegeneracy of the mapping e'^'=A, we get that the dimension of subspace 

{AJ{f)'ip, -^J{f)'ip ,..., -^J{f)'il)) 

CLTi CLTk 

equals the dimension M of the linear span of {k + 1) following vectors 

AJ{fc)'ip A ao{f)Ae~'"*‘'^gi,AJ{fc)i) ak{f)e~'"*='^gi, 

■^J{fc)'il) + ai{f)e~^'=^gi,-T^J{fc)'il) + ak-i{f)e~'"^^gi. 
dTi CLTk-l 

By inductive assumption there exists point fc such that the vectors 

AT(rc)V', -^Jijc)i2, ..., -T^Jijcl'il) 

CLTi dTk-l 

are linearly independent. By lemmas [5] and [6] there exists Tk ^ 0 such that the vectors 
AJ{fel'll), -^Jifel'll ),..., J{fel'll), Ae~'^^^gi,e~'^^'^gi 

dTi dTk-i 

are linearly independent and ao(f) = —2{e'^'^^J{fcl'ip,gil2 7 ^ 0. Thus ^ = k + 1 and induction 
step is proved. ▼ 






Lemma 5 For any vector 'ip ^ L \ {0} the set 


{t>0: gi)2 7 ^ 0} C M+ 

is open end everywhere dense subset. 

Proof. To prove this note the following equalities 

The second because A* = —A, where A* be the adjoint operator to A for the scalar product 
(, )h, that follows from the equality 

{Au,w)h = {Vp,q ')2 - {Vq,p ')2 = -{{Vq,p ')2 - (p, Vq') 2 ) = -{u,Aw)h. 
for u = {q,p)^, w = {qPp'Y ^ L. Then we use 

CO 

{'ip,e-^^gi)H = ^{-l)^ — ^,A^gi)H 

k=0 

as, by Lemma [T^ the linear span of vectors vectors A^gi, /c = 0,1, 2,... coincides with L. Then 
by analiticity of the left-hand part of the latter formula we get the proof. T 

Lemma 6 For any linearly independent vectors wi,... ,Wk & L, k < 2N — 1 the set 

T{wi ,..., Wk) = {t > 0 : dim{e~^^gi, Ae~^^gi,wi ,..., Wk) = k + 2} cR+ 

contains an open and everywhere dense subset. 

Proof. If fc < 2N — 2, then choose vectors Wk+i, ■ ■ ■, W 2 N -2 so that dim(t(;i,..., W 2 N- 2 ) = 
2N — 2. Note that 

T{wi ,..., W 2 N- 2 ) C T{wi, ...,Wk) 

That is why it is sufficient to prove the assertion of the lemma for the case k = 2N — 2. Thus 
hrther on we assume that k = 2N — 2. 

Consider the following real function on L; 

W{u) = dei[u, Au,wi,... ,W 2 n- 2 ]- 

Denote X = {u G T : W{u) = 0} the set of zeros of the function W. It is clear that 

T{w,,..., W 2 N- 2 ) = {t > 0 : lT(e-*V) 7 ^ 0} = {t > 0 : ^ X} 

It is clear that IT is a quadratic function. Thus X is a quadric. Further we dehne the 
canonical type of the quadric X. Denote 

Li = {Wi, . . . , W2N-2)- 

Consider also the orthogonal complement Lj; to Li (in the scalar product {,)h)- By dehnition 
of W, we have Li C X. Thus in some coordinates {ui,U 2 , ■ ■ ■ ,U 2 n) on L the function W will 
look like 

W{u) = aiu\ -|- 02 ^ 2 ) 


for some ai, 02 G M (as the basis for such coordinates one could choose the vectors wi,..., W 2 N- 2 , ui,U 2 , 
where Ui,U 2 are the approriate coordinates of the orthogonal compliment Lj-. As Li C X, then 
the form will not depend on the first 2N — 2 coordinates, and one can choose two coordinates 
in Lj- by method of Lagrange). 

If we show that there is exist a vector u such that W{u) 7 ^ 0, than it will follows that Oi and 
02 cannot be simultaneously zero. And consequently the canonical type of X should be one of 
three types: a 2N — 2-dimensional hyperplane, 2N — 1-dimensional hyperplane, the union of 
two 2N — 1 hyperplanes. Then applying the same argument as in lemma [5] we get the proof of 
lemma O 

Let show that that there is exist a vector u such that W{u) 7 ^ 0. 

For vector u ^ L denote {Au)^ the orthogonal projections on of the vectors u, Au 
correspondingly. As the determinant is polilinear we have 

W{u) = det[M-^, {Au)^,wi,... ,W 2 n- 2 ]- 

That is why X is the set of all u G L such that M-*-and [Au]^ are linearly dependent. Note that 
for any u = {q,p)'^ G L 

{u,Au)h = {Vq,p)2 - {p, Vq)2 = 0 . 

Consider two cases: 


1. Subspace Li is invariant with respect to A. Let us show then is also invariant with 
respect to A. In fact, for u G n G Li we have 

{Au,v)h = {u,A*v)h = -{u,Av)h = 0 , 

as Av G Li by invariance of Li. Then is invariant with respect to A and X fl is 
the set of all u G such that u and Au are linearly dependent. As the spectrum of A is 
pure imaginary and does not contain zero, then for all n 7 ^ 0 G the vectors u and Au 
are linearly independent. Then 

XAL^ = { 0 } 

and the proof is finished in this case. 

2. Subspace Li is not invariant with respect to A. In this case there exists vector n G Li 
such that u = (An)-*- 7 ^ 0. Consider vector u G such that u, u are linearly independent 
and u\u are linearly dependent, where u' G is a vector orthogonal to u in L^*-, again 
with respect to scalar product (, )h- As (m, Au)h = 0 , then for the projection (An)-*- of 
the vector Au on the subspace we have 

(Am)-*- = cu\ 


for some c G M. If c 7 ^ 0, then it is clear that W{u) 7 ^ 0. Assume that c = 0 and consider 
the vector 

u* = V + u. 


Then 

(m*)-*- = -u, (Am*)-*- = (An -f Ah)-*- = u + (Ah)-*- = u. 

By definition of h the vectors h and u are linearly independent. Then (n*)-*- and (An*)-*- 
are also linearly independent, and then W{u*) 7 ^ 0. The proof in this case is also finished 

▼ 


3.4 Proof of theorem 1 


Fix two points V'l) '02 ^ We want to show that for some m*, not depending on 0i, 02, there 
exist Ti,..., Tm* ^ 0 such that 

J(ri,... ,rm*)'0i = "02 

By local covering theorem there exists point f = (ti, ..., r^) G such that there exist 
neighborhood 0(f) C of f and neighborhood O('0i) C At of -0* = = J(ri,..., rm)0i 

such that 

J*(0(t)) = 0(V>:). (22) 

Consider e-neighborhood 0^(0]') of 05^ in the norm, corresponding to the scalar product (, )h, 
such that 

Oei'ipV) = {0 e At : 110 - 0 *||h < e} C O(0i) 

By closure theorem there exist 0,..., r0 ^ 0 such that 

-02|k < e. (23) 

Note that for hxed t[, ... the map J{t[, ..., r0) : L ^ L is non-degenerate and conserves 
the norm || \\h. That is why from (|23|) the following inequality follows 

1101 - J"^(r0.. .,rj 02 ||H < e. 

Thus, J“^(r(,..., r0)02 G 000*) and by (l2^ there is f* = ..., t^) G 0(f) such that 

Am)02- 


This can be rewritten as 


j • • • A0• • •, T-0)01 = j(0,..., r0) J(ri*,..., r 0 ) 0 i = 02- 


▼ 


4 Stochastic part - proof of convergence theorem 

4.1 Embedded process 

Here we consider the sequence ipk = which is a discrete time Markov chain (embedded 

chain) with state space Ai. We prove that it is (Markov) ergodic as it is dehned in the following 
theorem concerning more general class of discrete time Markov chains. Namely, we consider 
Markov chains on compact state space X with Borel cr-algebra B{X) (with countable basis) 
and transition probability kernels P{x,A), which are probability measures on B{X) for any 
X G X and measurable functions on X for any A (Z X. We will consider the class of such chains 
satisfying the following assumptions; 

Al) for some integer m > 1 and any x G X the m-step transition probability P"^{x,.) is 
equivalent to some hnite non-negative measure p such that /i(0) > 0 for any open set O Z X. 
Moreover, for any x there exists m-step transition density p'^{x,y) (with respect to /i), which 
is measurable on Ad x Ad; 

A2) for any open O Z X the function P{x, O) is lower semi-continuous. 

Theorem 6 Under assumptions Al, A2 the Markov chain is ergodic, that is there exists 
probability measure it on X such that 

sup \P^{x,A)-7r{A)\^0 (24) 

A&B{X) 


as n ^ oo, uniformly in x. 

This theorem follows probably from the existing deep theory of such Markov chains with 
continuous state space, see for example mm, but we did not hnd the statement we need, 
and for the reader’s convenience we give a short proof below, using the ideas from mm- 

Corollary 2 The embedded chain tjjk has the unique invariant measure 7i, and moreover it is 
ergodic as in theoreml^ 

We have only to prove the properties Al) and A2) for our embedded chain. For the rest of 
this section the integer m is the same as in the covering theorem. 

Proof of Al) We will prove that the measures tt and P™'('^, A) are equivalent for any fj. 

Lemma 7 For any measurable B Z Ai its Liouville measure vr(i?) = 0 iff the Lebesgue measure 
A of the set in Qrn is zero. 

Proof. 1) Assume that for some P C Ad we have vr(P) = 0. Let us show that X{{Jff)~^{B)) = 
0. Let Acr be the set of critical points of the map Jff (that is points r = (ri,..., r^) where 
the rank is not maximal) and let E = Jtf[Acr) C Ad be the set of critical values of 
By Sard’s theorem vr(P) = 0. But as = Ad, then there exists non-critical point 

r = (ti, ..., Tm) G kirn, that IS such that the rank of djff at this point equals 2N — 1. As the 
map Jfr is analytic in the variables ri,..., r^, the set of points where the rank is less than 
2N — 1, has Lebesgue measure zero. Then the equality \{{Jff)~^{B)) = 0 follows from theorem 

1 of m 

2) Assume that for some P C Ad we have 7 r(P) > 0, and let us show that A(( J,^)“^(P)) > 0. 
By Lebesgue differentiation theorem there exists point if' & A4 \ E and its neighbourhood 


0{'ijj') such that 7r(O('0') fl i?) > 0. Then there is point r = T{'ijj') G some 

its neighborhood 0{t) C Qm, so that the restriction of on 0(r) is a submersion. Then 
7^{0{^jJ') n B) > 0 implies fl 0(r)) > 0. ▼ 

Denote the product of m densities p, then as for any B G Ai 

P’"(i>,B)= ( p<”‘>(T)dr. 

by lemma [7] we get that P™ and vr are equivalent measures. 

The proof of measurability of the transition density there is in theorem 1, p. 180 of |16| . 
and in Proposition 1.1, p.5, of HU. 

Proof of A2) 

Lemma 8 Our Markov chain 'ipk O a weak Feller chain, that is for any open O C M the 
transition probability P{'if, O) is lower semicontinuous in fj- 

Proof. For any ip denote lp(r) the indicator function on P+, that is lp(r) = 1 if G O, 

and zero otherwise. Then we have 

P{'iP,0)= [ lp(r)p(r)dr, 

J M-|- 

Let xpn ipi'ipn ^ Ai. as n —)■ oo. Fix r ^ 0 and consider two cases: 

1. G O, then starting from some n the inclusion le'^^ipn ^ O holds, as O is open. 


That is why 

lim l^„(r) 

= ipW = 1 


n—>-oo 

2. ^ 0. Then 

lim inf (r 

n 

) ^ lp(r) = 0. 

Thus for any r 

lim inf 

n 

(r) ^ l^(r) 

Then by Fatou lemma 

lim inf Pi'ipn, 

n 

0)^P(V^,0) 


▼ 


4.2 Proof of theorem [6] 

Small sets We will need the following important dehnition. 

Definition 7 Below u will be any non-zero non-negative measure not necessarily probabilistic. 
The Borel subset C G X is called -small (or simply small, if it is -small for some 
integer n > 0 and some v ) if for any x & C and any Borel set B 

P^{x,B) > iy{B) 

Lemma 9 Assume 41). Then for some rii ^ 1 and some v there exists {u,ni)-small subset 
C G B{X) such that n{C) > 0. 


Proof. It follows from the assumption Al) that for any n ^ m there is measurable function 
Pn(x,y) such that 

P^{x,B)= [ pnix,y)p(dy), 

J B 

for any x E X, Be B{X). Moreover, Pni^x, y) > 0 for almost any {x,y) E X x X (with respect 
to /i X /i). 

We will prove more. Namely, that for some n ^ m there exist sets Bi, B 2 E B{X) of positive 
measure p, and some constant 5 > 0 such that for all a; G i?i, y E B 2 we have 

Pn{x,y)>5 (25) 

As the density pm{x,y) is measurable and almost everywhere positive, one can find number 
c > 0 so that the sets 

A = {{x,y) E X X X : Pm{x, y) > c] 

{{x^y^z)EXxXxX\ [x^y) E A VL (?/, z) E A} = [Ax X) fi (X x A) 

have positive measures p^ = pxponXxX and p^ = pxpxponXxXxX correspondingly. 
Denote Or(x) C X the open neighborhood of a; G X of radius r, and put Or{x,y) = Or{x) x 
Or{y) C X X X. By Lebesgue differentiation theorem there exists a set Aq of zero measure 
p X p such that for any (a;, y) G A \ Aq 

p^{AnOr{x,y)) ^ 
p\Or{x,y)) 

It follows, that the set ((A \ Aq) x X) fl (X x (A \ Aq)) also has positive measure. Consider 
some point {x*,y*,z*) in this set. Choose r so that the following inequalities hold; 

p\AnOr{x*,y*)) > ^p\Or{x*,y*)) 
p\AnOr{y*,z*)) > ^p\Or{y%z*)) 

For any x E X put 

^dx) = {y E X ■. {x,y) E A}, Ar{z) = {y E X ■. {y, z) E A} 

Define Bi by 

= {a; G 0,.(a;*) : p{Al{x) AOr{y*)) > ^p{,Or{y*))} 

Otherwise speaking, the set Bi consists of the points x E Or{x*), for which the set Al{x) is 
sufficiently large inside Or{y*)- Let us show that Bi has positive measure. Assume the contrary 
- that for almost all points of Or{x*) the inequality p{AL{x)AOr{y*)) ^ \dOr{y*)) holds. Then 
by Fubini theorem 

p^{AnOr{x\y*)) = [ p{AL{x)nOr{y*))p{dx) ^ ^p{Or{y*))p{Or{x*)) = ^{pxp){Or{x*,y* 
JOr{x*) 4 4 

That contradicts the choice of the points x*,y*. Thus, p{Bi) > 0. Similarly, one can show that 
the set 

B 2 = {zE Or(z*) : p{Ar{z) n Or{y*)) > -MOr{y*))} 



has positive measure /r. But from the dehnition of the sets -Bi, B 2 it follows that for any points 
X ^ Bi, z ^ B 2 the following inequality holds 


n Ar{z)) ^ 

Also we have the estimates for the density p 2 m and any x & Bi, z ^ B 2 

P 2 m{x,z)= / Pm{x,y)pm{y,z)p{dy) ^ / Pm{x,y)pm{y,z)p,{dy) > 

Jx J Ai^{x)r\AR{z) 

> p{Al{x) A Ar{z)) > ‘^p{Or{y*)). 

Thus we have proved (|25|) . Now we finish the proof of lemma |9l Note that for any x G -Bi 
and any B G B{X) we have 

P”(x,5) = f pn{x,y)p{dy) ^ [ Pn{x,y)jj,{dy) ^ 6 p{BnB 2 ). 

Jb J BrB2 

Moreover, as the sets Pi, B 2 have positive measure, then there exists subset C <Z B 2 oi positive 
measure p and constant 5' > 0 such that P‘^{x,Bi) > 6 ' for all x E C. It follows that for all 
X G C, Be B{X) the following inequalities hold; 

F™+"(x,P)= [ P^ix,dy)P^iy,B)^ [ P^ix, dy)P^iy, B) ^ 6 ' 6 piB n B 2 ). 

Jx J Bi 

Thus, C is {i',n + m)-small, where ^{B) = 66 'fi{B fl B 2 ), and moreover i^(C') > 0. ▼ 

Lemma 10 Let C E B{X) be (z/, ni)-small and assume that for any x from some set D E B{X) 

P'^{x,C) > 6 . 

for some 5 > 0 and n ^ 1. Then the set D is { 6 u, n + ni)-small. 

Proof. For any x E D and any B E B{X) by semigroup property 

P^^^^{x,B)= [ P^{x,dy)P^^{y,B)^ [ P^{x,dy)P^^{y, B) ^ P^{x,C)u{B) ^ 6 u{B), 

Jx Jc 

▼ 

Lemma 11 Under assumptions Al) and A2) the set X itself is small. 

Proof. Consider some (z/, ni)-small subset C E B{X) of lemma |9l For k = 1,2,... introduce 
the subsets 

Ak = {xEX: P^{x,C)>j}, 

k 

where m is dehned in Al). From lemma fTOl we have that A^ is a + ni)-small set. Let 

us prove that its closure is also a small set. Note that the measure v is regular, that is for any 
Borel set B 


z/(P) = sup{z/(A:)}, 


where the supremum is over all compact sets K <Z B. As v{C) > 0, there exists compact set 
K <Z C such that for all x ^ 

^ = 5 > 0 . 

k 

Let {xn\n=i, 2 ... £ Ak and —)■ x, then by semi-continuity of the transition probability we get: 

= l-P^+^^{x,X\K) ^ l-\imm{P'^+^^{xn,X\K) = limsup P”"+”i(x„, iL) ^ 

Thus, for any x G we have P^~^^^{x,K) >6/2, and applying lemma fTOl we get that A^ is 
a (I z/, 712 )-small set for some n 2 . 

But also by assumption Al 

^ = ur=iA, 

As X is not a countable union of sets which are nowhere dense, then for some k there is an 
open subset O C A^. Thus for any x E X and any B G B{X), we have 

pm+n2^x,B) = J P^{x,dy)P^^{y,A)^ j P^{x,dy)P'^\y,B)^^-v{B)P^{x,0). (26) 

As the measures P(x, •) and y are equivalent for any x G X, we have P™'(x, O) > 0. But 
P™(x, O) is lower semi-continuous, and thus attains minimum on the compact X. It follows 
that 

P"‘(x,0) > P 

for some P > 0 and any x E X. Using inequality fl26l) . we get the proof. 

Proof of theorem |6] ForA G B{X) and n ^ 1 denote 

4(A) = inf P^x, A), 4(A) = supPAx, A). 


Note that 

4+1 (A) = inf f P{x,dy)P^{y, A) ^ ini f P{x,dy)In{A) = In{A). 

xGX xGX 

Thus,, for hxed A G B{X) the sequence 4(A) is non-decreasing. Similarly the sequence 4(A) 
is non-increasing. We shall prove that 4(A) — 4(A) tends to zero as n —)■ cxd. 

Take number N and measure v as in Lemma [11] Then for any rz, ^ 1 

P"+^(x,A)= [ P^{x,dy)P-{y,A)= [ {P^{x,dy) - u{dy))P-{y,A) + [ u{dy)P-{y,A) 

Jx Jx Jx 

As the measure P^(x, •) — z/(-) is non-negative, we get from this equality that 

P"+^(x,A)^4(A) [ iP^ix,dy)-u{dy))+ [ z/(4)PA7/, A) = (l-z/(X))4(A)+/ uidy)P-{y,A) 
Jx Jx Jx 

and then 

4+iv(A) ^ (1 - h)4(A) -f Cl, <5 = z/(X), Cl = / iy{dy)P^{y, A). 

Jx 


Similarly we get the upper bound for Sn+N{,A): 

A) < (1 - iy{X))Sn{A) + I v{dy)P^{y, A) < (1 - + ci 

and 

Then the difference is estimated as follows 


5„+;v(dl) - In+N{A) ^ (1 - <5)(^„(A) - 4(A)). 

From the last inequality and monotonicity of the corresponding sequences it follows that for 
any x E X there exist the following limits and that they are equal 

lim 4(A) = lim Sn{A) = lim Pn{x,A) = vr(A), 

n^oo n—>-oo n—>-co 

and moreover the convergence is uniform in A G B{X) and in x. ▼ 


4.3 Proof of theorem [2] 

We will use the following theorem (strong law of large numbers) for discrete time Markov chains 
on arbitrary state space X equipped with cr-algebra A. Let P'^{x,B) be n-step transition 
probability assumed to be measurable on X for any B E A and is a probability measure on 
{X, A) for any x. Let us assume that there exists invariant measure tt on (X, A) such that 
uniformly in x 

sup |T’"'(x, A) — 7r(A) | —)■ 0, n —>■ cx) 

A&A 

Denote Px the measure on trajectories (xq = x,xi,X 2 ,...) with initial point x. Under these 
conditions the following assertion holds. 

Theorem 8 For any f E L^(X, tt) and any x E X we have Px-a.s. 

1 A f 

hm-y'/(xfc)= / f{x)7r{dx) 

Proof. See [H], p. 140, and [l6], p. 209. 

To prove theorem [2] we need the following lemma. 

Lemma 12 For any measurable bounded function f on M. and any initial state '0(0) = "0 the 
following limit holds a.s. 

lim Mfitn) = vr(/) 

N^oo 

Proof. Denote X^ = (0^, t^+i), fc = 0,1,... the Markov chain with values in X = AT x M_|_. 
Then 

/ f{m)ds= / f{e^-^-^^^Uk)ds= / f{P^^l^k)ds = F{Xk) (27) 

J J Jo 


where 


Then 


1 ^ rtk+i 1 ^ ^ 

= -J2 /(V’(s))* = -J2 (28) 

i.o t-o 

It is easy to show that Xj. has invariant measure /i = tt xP^, = pds, satisfies the conditions 

of theorem [S] as 'ipk satisfies it. Then 

Af-l 


where 


N^oo iV ^ / V 

fc =0 


p{F) = I Pr{dt) I 7r{di/j) I dsf{e^^'tp) = I Pr{dt) I ds I 7r((i-^)/(e^"^V’) = 




IM 


'0 J M 


= 7r(/) / Pr{dt) ds = 7r(/)Pri 

J M_|_ J 0 

Moreover, by strong law of large numbers for independent random varoables we have 

lim ^ = Eti 

N^OO N 

Then by fl28l) we get the proof of the lemma. ▼ 

To prove theorem [2] we have to estimate the difference between Mf{t) and Mf{t]s[). Using 
the boundedness |/('?/’) | ^ c we have 


\Mf{t) - Mf{tN)\ ^ It / f{'ip{s))ds\ + 


|t — tpf\ 


t 


tjv 


t 




|t — 


[c+\MfM\). (29) 


For any t > 0 define the random index N{t) so that 

tN{t) ^t< t]\fp)+l- 

and note that 

It — tNipl . TN{t)+l _ TN{t)+l 


* ' *'V„| zSn 

As Eti < oo, the law of large numbers, as iV —)■ oo, gives a.s. 

T'Ar+l 


Z^fc=i p 

But N(t) —)■ cx) as t ^ oo. Then the right-hand side of (|2^ tends to 0 a.s. as iV = N(t) and 
t —)■ oo. ▼ 


^ 0 








5 Appendix 

5.1 Proof of lemma [2] 

For any t ^ 0 one can show that 

tA ^ f cos{t\/V) {Wy^ sm{t\/V) \ 

\ —\/V sin(tVl7) cos{t^/V) ) ’ 


where \/V is the positive square root of the matrix V. Then for any k 

= cos{ukt)Qk - UkSm{ukt)Pk, 

e^Apj, = + cos{ukt)Pk, 

OJk 


1 ,..., and f ^ 0. 


(30) 

(31) 


Using iOD-dslD we have for any 

N , 

ipit) = ^ icos{ukt)qk + 

k=l ^ 


sm{ukt) _ \ 

- Pk 

OJk J 


N 

Qk + '^ {-OJk sm{ukt)qk + cos{ukt)pk) Pk- 

k=l 


and then 


~ N / X ~ sinfoifcf) _ 

qk{t) = cos{Ukt)qk H- Pk 

OJk 

Pkit) = - a;fcsin(a;A:t)gfc + cos(a;fcf)pfc. 

From these two formulas we see that pair of functions {qk(t),pk(t)), k = 1,... ,N, corresponds 
to the dynamics of one-dimensional oscillator of unit mass and frequency oJk, where qk{t) is the 
oscillator coordinate and pk(t) is its momentum. Thus the dynamics is isomorphic to the 
uniform movement on the torus with velocity {uf,... jOjjf). This gives the second assertion. 

5.2 Case N = 1 

Note the identity for any pj 

= e-^^^Ipj (32) 

Let us hrst prove the proposition for t = Tq, p) = pjQ = (1,0) and arbitrary p)'. Then for any 
dehne the one-to-one mapping W : (0, |To) —)■ S. namely for any ti G (0, \Tq) we dehne 

W{t') =P)' = 

Then it is sufficient to take t' = and choose minimal t" > 0 so that 

pj’ = e^^"pj. 

For t > To the proof is quite similar but we will not need this case to prove convergence. 






5.3 Mixing subspace 

Here we give some properties of the mixing space. The following two lemmas show how the 
dimension of L_can be explicitely characterized. 

Lemma 13 The space L_ is invariant with respect to A. Moreover 

L_ = ({HV: fc = 0,l,...}) 

where ( ) is the linear span of the set of vectors. 

Consider also the orthogonal complement to L_ in the scalar product (, )2 

Lo = L\ 

Then it is also invariant with respect to A. Moreover, the vector -0 G Tq iff for the hamiltonian 
dynamics with initial condition = 'if{0) the momentum pi{t) = 0 for any t. 

Lemma 14 Assume that the spectrum ofV is simple, and let {ui,... ,V]s[} be the eigenvectors 
of V, they form a basis in . Then the dimension of Lq is twice the number of Vk having 
coordinates Vk,i = {ei,Vk) = 0. 

What occurs if condition (jl]) is not fulhlled, exact formulas for the dimension of L_ for the 
chain of harmonic oscillators and for other cases see in [IIl[inil9]. 


5.4 Proof of lemma E] 

We have 


\v*k 





Ay 


\Pk 


+ 27n V7n + c(l - 7n) + c(7fc " ll) 
Ay 


As for any A; = 1,..., we have '^k ^ '^n ^ 0, then the expression under module in the last 
formula is non-negative, and we have 


\Pk\ = l/^fcl [y-c 


Ay 

Consider two cases: 

1. 7 fc ^ 1. Show that f~^{x) is monotone increasing, that is its derivative 


(/^(^))'=o 1 + 


a;(l — c) 


-f c(l — x‘^) 


> 0 . 


Thus y = f^ipin) ^ f^irtk)- Taking into account ^ 1, we have the inequalities 


\Pk\ = l/^fcl (y-c 


= IA 


.1 -7^ 

Ay 


^ l/^fcl f Aik) 


| 27fc + 27fcV7fc + c(l-7fcy 

4/+(7fc) 


c(1-7D \ 

4/+(7fc) / 


\ldk\lk = Tk- 
















2. 7 fc > 1- Then we will show that / ( 7 ^) ^ y. It is easy to check that / (x) is increasing 
and then 

0 < /“(7fc) ^ /“(Ttv). 

Note that the left inequality holds because /~(1) = 0. 

Let us prove that for any x G M we have f~{x) < |x. Consider the chain of inequalities 


/ (^) - ^ (^(1 - c)x - 7 /x 2(1 -c) + c) 


1 / -c - x^jl - c)c \ ^ ^ 

2 y (1 - c)x + ^Jx‘^{l - c) + cJ 


Thus we have proved that 

/■(7fc) ^ /“(7v) < ^7iv. (33) 

By /"'"(O) = |\/c and the evident inequality (/+(x))' ^ which holds as c ^ 1, we have 


y = > ^(c7n + \/c)- 

Then by ((331), (El]) and c ^ 


(34) 


y - / (7fc) ^ ^(-c(7w - 7n) + \/c) = ^ (-\/c(7iv - In) + 1 ) ^ 0. 

Remind that we have proved earlier that the function f~^{x) is increasing, then taking 
into account the latter inequality we have the following bounds 


Kl = lAI ( V - cLji) = l&l {v + cT_!l) ^ 1 , 3,1 7 +( 7 ,) + c-T 1 


= l& 


Ay J V 

Ik-Ik- c(l - 'll) + c(7l - 1) 


4/“(7fc) 
= |/5fc|/’^(7fe) ^ l/^fc|7fc = rk- 




4/ (7fc). 

2(71 - l)/+(7.) 
2c(7| - 1 ) 


In the last inequality we used fact that f~^{x) < x for x > 1 holds due to /+ + / = x. ▼ 
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